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Grand unified theories can lead to non-universal boundary conditions for the gaugino masses at 

O ! 

the unification scale. We consider the implications of such non-universal boundary conditions for 
Cs| ' the composition of the lightest neutralino as well as for the upper bound on its mass in the simplest 

§: ^ „ 

neutralino and chargino masses in different representations of SU(5) which lead to different non- 
universal boundary conditions for the gaugino masses at the unification scale. We also consider the 
phenomenological implications of the non-universal gaugino masses arising from a grand unified 
theory in the context of Large Hadron Collider. In particular we investigate the detection of heavy 
neutral Higgs bosons H°, A from H°, A — ► X2X2 ~^ 4Z, and study the possibilities of detecting 
the neutral Higgs bosons in cascade decays, including the decays X2 ~^ h°(H°, A )x± — ► bb\\- 
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I. INTRODUCTION 



Supersymmetry is at present an attractive framework in which the Higgs sector of the 
Standard Model (SM), so crucial for its consistency, is technically natural. It is widely 
expected that some of the supersymmetric partners of the SM particles will be produced at 
the CERN Large Hadron Collider (LHC) which is going to start operation in a few years time. 
In the experimental search for supersymmetry (SUSY) the lightest supersymmetric particle 
will play a crucial role since the heavier supersymmetric particles will decay into it. In 
SUSY models with R-parity conservation, the lightest supersymmetric particle is absolutely 
stable. The lightest supersymmetric particle is constrained to be a weakly interacting neutral 
particle 

In most of the supersymmetric models the lightest neutralino (x§), which is typically an 
admixture of gauginos and higgsinos, is the lightest supersymmetric particle (LSP). Such 
an LSP is a good candidate for a particle dark matter |2j. From the point of view of 
experimental discovery of supersymmetry at a collider like the LHC, the LSP is the final 
product of the cascade decay of a SUSY particle. In this work we will assume that the 
LSP is the lightest neutralino, and that it escapes the collider experiments undetected. 
The cascade chain will typically also contain other neutralinos (x®, j = 2,3,4) as well as 
charginos (xf , i = 1,2). The charginos are an admixture of charged gauginos and charged 
higgsinos. The composition and mass of the neutralinos and charginos will play a key role in 
the search for supersymmetric particles. These properties determine also the time-scale of 
their decays. The mass patterns of the neutralinos in models with different particle content, 
or with specific SUSY breaking patterns were considered in some detail in 3, 0| . 

Although most of the phenomenological studies involving neutralinos and charginos have 
been performed with universal gaugino masses at the grand unification scale, there is no 
compelling theoretical reason for such a choice. Gaugino masses follow from higher dimen- 
sional interaction terms which involve gauginos and auxiliary parts of chiral superfields in 
a given supersymmetric model. Assuming an SU(5) grand unified theory (GUT) model, 
the auxiliary part of a chiral superfield in these higher dimensional terms can be in the 
representation 1, 24, 75, or 200, or some combination of these, of the underlying SU(5) 
gauge group. If the auxiliary field of one of the SU(5) nonsinglet chiral superfields obtains 
a vacuum expectation value (VEV), then the gaugino masses are not universal at the grand 
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unification scale. Moreover, nonuniversal soft supersymmetry breaking masses, like gaugino 
masses, are a necessary feature in some of the supersymmetric models, e.g. in anomaly 
mediated supersymmetry breaking models the gaugino masses are not unified Q]. 

As indicated above, the phenomenology of supersymmetric models depends crucially on 
the composition of neutralinos and charginos. Thus, it is important to investigate the 
changes in the experimental signals for supersymmetry with the changes in the composition 
of neutralinos and charginos that may arise because of the changes in the underlying bound- 
ary conditions at the grand unification scale, or when the underlying supersymmetric model 
is changed. The implications of nonuniversal gaugino masses has been considered in a num- 
ber of works, e.g. in a study of constraints arising from exp eriment al measurements 0,0,0], 
and in the context of supersymmetric dark matter 

lightest neutralino were studied in the context of nonuniversal gaugino masses. 

In this paper we shall study the implications of the nonuniversal gaugino masses for the 
phenomenology of neutral Higgs bosons. It has been known for quite some time that the 
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cascade decays of the SUSY particles may be a major source of the Higgs bosons 
the copiously produced strongly interacting particles can cascade decay to the Higgs bosons. 
In addition to the obvious interest in producing the Higgs bosons, it has been realized that 
this method of producing the Higgs bosons does not depend on the value of tan (3. Thus, this 
method of producing Higgs bosons may help to cover a larger parameter space as compared 
to the more conventional methods of studying the Higgs sector of supersymmetric models, 
including also the heavier Higgs bosons. The gauginos also play an important role in the 
decays of Higgs bosons when they are kinematically allowed to decay to the second lightest 



inos and two leptons 
I5I 16^ . We note here that 



neutralino pair, which in turn may decay to the lightest neutra 
Such a signal seems to be relatively easy to discover at the LHC 
Higgs boson production via cascade decays and detection via Higgs decay to neutralinos 
has been studied in CMS detector simulations at LHC Q, 3] in the case of minimal 
supersymmetric standard model (MSSM) with universal gaugino masses. Here we study the 
Higgs production and decay when gaugino masses are nonuniversal. 

The plan of the paper is as follows. In Section II we consider in detail the nonuniver- 
sality of gaugino masses as it arises in SU(5) supersymmetric grand unified theory. In this 
Section we consider analytically the implications of such a nonuniversality for neutralino 
and chargino masses. We derive sum rules involving the neutralino and chargino squared 
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masses when the supersymmetry breaking gaugino masses are nonuniversal. In Section III 
we consider the phenomenology of Higgs bosons when the gaugino masses are nonuniversal. 
In this Section we consider Higgs decays to heavier neutralinos which then cascade into the 
lightest neutralino and leptons. In Section IV we calculate the production of squark and 
gluino pairs in a particular scenario where the gluinos are heavier than squarks, and then 
study the cascade decays of the squarks into Higgs bosons. We conclude our paper with a 
summary in Section V. 

II. NONUNIVERSAL GAUGINO MASSES IN SUPERSYMMETRIC SU(S) 

The masses and the compositions of neutralinos and charginos are determined by the soft 
supersymmetry breaking gaugino masses Mi, M2, and M3, corresponding to U(l), SU(2), 
and SU (3) gauge groups, respectively, the supersymmetric Higgs mixing parameter /x, and 
the ratio of the vacuum expectation values of the two neutral Higgs bosons H® and H®, 
(H®) I (Hi) = tan (3. In the simplest supersymmetric models with universal gaugino masses, 
Mi, Mi, and M 3 are taken to be equal at the grand unified scale. However, in supersymmetric 
theories with an underlying grand unified gauge group, the gaugino masses need not be equal 
at the GUT scale. In this Section we consider the nonuniversality of gaugino masses as it 
arises in the simplest of the supersymmetric grand unified theories, namely supersymmetric 
SU (5) grand unified theory, and its implications. 

In grand unified supersymmetric models, including SU(5) grand unified models, non- 
universal gaugino masses are generated by a non-singlet chiral superfield $ n that appears 
linearly in the gauge kinetic function /($) (the chiral superfields $ are classified into a set 
of gauge singlet superfields $ s , and gauge nonsinglet superfields <3> n , respectively under the 
gra nd unified group), which is an analytic function of the chiral superfields $ in the theory 
|l8| . If the auxiliary part F$ of a chiral superfield <3> in /(<&) gets a VEV, then gaugino masses 
arise from the coupling of /($) with the field strength superfield W a . The Lagrangian for 
the coupling of gauge kinetic function to the gauge field strength is written as 



where a and b are gauge group indices, and repeated indices are summed over. The gauge 




(1) 
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kinetic function f a b(&) is 

frtW = fomSab + E /n(^)F + ■ • ( 2 ) 

n M P 

where as indicated above the $ s and the $ n are the singlet and the non-singlet chiral 
superfields, respectively. Here fo($ s ) and f n (& s ) are functions of gauge singlet superfields 
$ s , and M P is some large scale. When Fq> gets a VEV (F$), the interaction ([I]) gives rise 
to gaugino masses: 

C a j.. D %^A a A b + /i.c., (3) 

where A a,fc are gaugino fields. Note that we denote by Ai, A2 and A3 the £7(1), SU(2) 
and 577(3) gauginos, respectively. Since the gauginos belong to the adjoint representation 
of 577(5), $ and F$ can belong to any of the following representations appearing in the 
symmetric product of the two 24 dimensional representations of SU(5): 

(24 <g> 24) Symm = 1 © 24 © 75 © 200. (4) 

In the minimal, and the simplest, case $ and F$ are assumed to be in the singlet repre- 
sentation of £77(5), which implies equal gaugino masses at the GUT scale. However, as is 
clear from the decomposition P|). $ can belong to any of the non-singlet representations 
24, 75, and 200 of SU(5), in which case these gaugino masses are unequal but related to 
one another via the representation invariants |19^]. In Table U we show the ratios of resulting 
gaugino masses at tree-level as they arise when F<j, belongs to various representations of 
577(5). For definiteness, we shall study the case of each representation separately, although 
an arbitrary combination of these is obviously also allowed. 

These results are consistent with the unification of gauge couplings 

qG = a G = a G = a G^ 1//25 ^ ( 5 ) 

at the grand unification scale, where we have neglected the contribution of nonuniversality 
to the gauge couplings. Such contributions have little effect on the phenomenological aspects 
that we are interested in this paper. Because of the renormalization group (RG) evolution 
we have at any scale (at the one-loop level) 



Mitt) = M t (GUT) 

ati(t) ai(GUT) ' 1 ' 
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TABLE I: Ratios of the gaugino masses at the GUT scale in the normalization M^{GUT) = 1, 
and at the electroweak scale in the normalization M^(EW) = 1 at the 1-loop level. 





Afp M 2 G M§ 


M EW 


M 2 EW 
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1 1 
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0.14 


0.29 
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24 


-0.5 -1.5 


1 


-0.07 


-0.43 


1 


75 


-5 3 


1 


-0.72 


0.87 


1 


200 


10 2 


1 


1.44 


0.58 


1 



Thus, at any scale we have 

5 a / Mi (GUT) \ q / M 2 (GUT) \ / M 3 (GUT) \ 

1 3cos 2 ^ V ai(GUT) / ' 2 sin 2 6 W \ a 2 (GUT) )' 3 3 ^ a 3 (GUT) j ' 

(7) 

For the 24 dimensional representation of 577(5), we then have 



M 1 1/5 a \ / 1 \ M 2 3 / a \ / 1 



M 3 " 2 V3 cos 2 ^iy/ Va 3 / ' M 3 2 Vsin 2 ^/ va 3 
Similarly, for the 75 dimensional representation of S77(5), we have 



(8) 



M 3 5 Gcos 2 ^) (q 3 ) ' M 3 3 (sin 2 ^) ("a)' (9) 

and for the 200 dimensional representation of 577(5) we have 

Mi_ /5 a \ / 1 \ M 2 / a Wl\ 

= Ucos 2 %J WJ ' M3 = Uin 2 WJ ' ^ ^ 

We can scale down these results to the the electroweak scale by using the relevant renor- 

malization group equations. At the electroweak scale we have the result Mi < M 2 for the 

singlet representation, \M\\ < |M 2 | for the 24 and 75 representation, and Mi > M 2 for 

200 dimensional representation of 577(5), respectively. The approximate values for the soft 

gaugino masses at the weak scale, Mi(EW) are shown in Table These are calculated using 

one loop RG-equations for the gaugino masses and the gauge couplings. Two-loop effect is 

to increase the Mi/M 2 -ratio slightly. 

In Fig.^we have shown the dominant component of the lightest neutralino (LSP) for the 

four representations as a function of tan/5 and M 2 (EW) for the value of soft supersymmetry 
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FIG. 1: Main component of the lightest neutralino in different representations of SU(5) that arise 
in the product @ for a common universal scalar mass = 1 TeV given at the GUT scale. The 
value of M2 is calculated and plotted at the electroweak scale. 

breaking scalar mass mo(GUT) = 1 TeV. The values of /1 used in the computations were 
determined by requiring the radiative electroweak symmetry breaking at the relevant scale. 
The universal trilinear coupling A was set to zero at the GUT scale and the sign of the 
Higgs mixing parameter /1 was set to +1, but the choice of the sign is not crucial to the 
composition of the lightest neutralino. The scan was done using the program SOFTSUSY 



2l\ that uses two loop RG /3-functions for the relevant parameters. 

For the case of the singlet representation, the dominant component is always the bino, as 
expected. This is also true for the 24 dimensional representation of SU(5). For the singlet 
case the experimental mass limit of the lighter chargino (m~± > 103 GeV if > 200 GeV, 
m~± > 45 GeV if mc < 200 GeV j^]) restricts the lower end of the M 2 range. In the 
24 dimensional representation the lower end of the M 2 range is restricted by the lightest 
neutralino mass limit m^o > 36 GeV j^j. 
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FIG. 2: The upper bound for (a) the lightest neutralino mass and (b) for the second lightest 
neutralino mass for different representations that arise in Q. 

For the 75 dimensional representation of SU(5), we have several possibilities. For the 
value of the soft parameter m = 1 TeV, one has a bino LSP for small values of M 2 , a 
wino LSP for slightly larger values of M 2 , and a higgsino LSP for M 2 ^ 300 GeV, all for 
a value of tan f3 £ 10. In the case of 75 dimensional representation there exists a band of 
discontinuity in the (M 2 , tan/3)-parameter space. For these values of parameters the lighter 
chargino mass becomes too light. The lower end of the M 2 range is restricted in this case 
by the experimental limit on the gluino mass. 

As seen in Fig. for the 200 dimensional representation the LSP is either a wino or 
a higgsino, depending on the values of M 2 and tan/3. Here, as in the singlet case, the 
experimental lightest chargino mass limit restricts the lower end of the M 2 range. Also in 
the 200 dimensional representation there is a small region around 7 £ tan (3 £ 8, 610 GeV < 
M 2 (EW) < 620 GeV, where the experimental mass limits of charginos (and also neutralinos) 
are not met. 

We recall here that there is a general upper bound on the mass of the lightest neu- 
tralino (x?) that follows from the structure of the neutralino mass matrix DO . This upper 
bound can be written as 



M|o < - [Ml + Ml + M| - sJ{Ml - Ml) 2 + M| - 2(Mf - Mf )Mf cos 26 w ) . (11) 

In Fig. |21 we plot this upper bound for the lightest neutralino mass for the four different 
representations of SU(5i) that we have considered in this paper. From Fig. |21 (a) we see 
that the large coefficients in the Table H] result in large differences in the upper bound on 
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FIG. 3: The sum rule 2J2 m.± — X) 

Tfl ~ o clS cL function of the gluino mass. 

Xi Xi 

the mass of the lightest neutralino for the four different representations in (j3J). Similarly, as 
n 

discussed in |4j, an upper bound can be obtained for the second lightest neutralino. This 
upper bound for \2 is shown in Fig. |2l (b). The gaugino masses here are calculated in the 



next-to-leading order (see e.g. "\). 

In order to study analytically the implications of the nonuniversal gaugino masses on the 
neutralino and chargino mass spectrum, we consider the trace of the neutralino and chargino 
mass squared matrices. From the trace of these matrices, we can calculate the average mass 
squared difference of the charginos and neutralinos. This mass squared difference depends 
only on the physical masses, and not on the Higgs(ino) mass parameter fi or the ratio of 



VEV's, tan/? 



201 ]. For the four different representations of SU(5) which arise in (£Q), we find 



at the tree-level the sum rules 



M 2 sum = 2(M 2 f + M|) - (M| + M| + M| + M|] 



M, 



2 



{a 2 2 -al)-^ + Am 2 w -2m 2 z , fori, (12) 



A at--al)^ + 4m 2 w -2m 2 z , for 24, (13) 



,9 o 1 o.Mj 

4 

M- 

{Va 2 2 -2ha 2 )—l r + Am 2 w -2m 2 z , for 75, (14) 

^3 

M~ 

{Aa 2 2 - 100a?)— f + 4m 2 w - 2m|, for 200. (15) 

a 3 
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From these sum rules we see that at the tree-level the average mass squared difference 
between charginos and neutralinos is positive for the representations 1, 24 and 75, whereas 
for the representation 200 it is negative. In this respect the representation 200 resembles the 
anomaly mediated supersymmetry breaking scenario, where it was found that the average 
mass squared difference is negative 



23[ . In Fig. we have plotted the above sum rules for 
the different SU(5) representations that arise in (0J). For the numerical evaluation of the 
masses, we have used the program SOFTSUSY |21[, including radiative corrections to the 
neutralino and chargino masses. 



III. HIGGS DETECTION USING H°, A -> X2X2 -» ^ 

It is often assumed, when considering the detection of the Higgs bosons in supersymmetric 
models, that supersymmetric partners are too heavy so that Higgs bosons cannot decay into 
supersymmetric particles. However, it may well be that for the heavy Higgs bosons H°, A , 
and H ± the decays to supersymmetric particles are important or even dominant On 
the other hand, the decay branching ratios of neutralinos and charginos have been analyzed 
in (24]]. In the case of large tan/5, when the couplings to the heavy fermions are enhanced, 
the decays to the third generation particles have been discussed in J^. For large values 
of tan/?, the decays to the third generation particles for the nonuniversal gaugino masses 
were discussed in |8|. Here we are interested in Higgs decay to which in turn decays to 
electrons and muons in the case of non-universal gaugino masses. 



A. Decay of xi *° leptons 

Of the supersymmetric particles, the light neutralinos X121 the light chargino Xii an d 
the lightest sleptons are usually among the lightest particles in the spectrum. Higgs decays 
to sleptons are suppressed because of the small coupling, which is proportional to the cor- 
responding lepton mass. The decay to the lightest neutralino LSP is among the invisible 
decays, which may be extremely difficult to detect at the LHC In the minimal SUGRA 
model, the second lightest neutralino and the lighter of the charginos have similar mass. In 
[15I ] , the decay of the heavy neutral Higgs boson to a pair of the second lightest neutralinos 
was studied. It was found that in the case when the branching ratio of x\ to two leptons 
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FIG. 4: Branching ratio as a function of tan (3, in the case of mj > m^,o and for the representations 
1 and 75. 

and the lightest neutralino is large, then the possibilities of detection are promising. Even 
though the branching ratio to a chargino pair may be larger jl^ . l^fil ]. the decay to x^ s is 
more promising because of the clear four lepton signal. We will, therefore, study the decay 
chain 



H°,A° 



X2X21 X2 



Xil + l , I = e,n 



(16) 



for the four different representations of SU(5) in (J3j). The decay X2 ~~ ^ depends on 

the parameters M 2 , M 1; /i, and tan/?, which control the neutralino masses and the mixing 
parameters, and also on the slepton masses m\. As long as the direct decay of \\ i n t° + 
is suppressed and the sleptons are heavier than the three body decays of X2 charged 
leptons and Xi wn l be significant. There can also be constructive or destructive interference 
between the Z° and the slepton exchange amplitudes which can have strong influence on the 
branching ratio. In some cases, we also consider the possibility that the sleptons mediating 
the decay X2 ~^ Xil + l~ can be on mass-shell. Also, in order to have large branching ratios to 
the sleptons, they should be lighter than the squarks. This is usually true in SUSY models. 



In the singlet case, the three body decay X2 ~^ Xil + l was discussed in Ref. jl5| for a 
particular set of parameters, for which the branching ratio is large. Here, all the sleptons 
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(including the staus) are assumed to have soft SUSY breaking masses of 250 GeV and 
the value of fi = -500 GeV. In this analysis M2 is a free parameter, and its value at the 
electroweak scale is taken to be 150 GeV. The squark masses are all taken equal to 1 TeV. 
We have also taken a large value of the trilinear scalar coupling A t — 1 TeV in order to have 
experimentally acceptable mass for the lightest Higgs boson. All the soft scalar masses, 
the value of A t and the value of /1 are taken at the electroweak scale. The pseudoscalar 
Higgs mass rriA is a free parameter and its value is taken to be 340 GeV. The value of Mi 
is determined from the ratio of the gaugino mass parameters in the singlet representation 
of SU(5) in (0J). Due to the mentioned choice of M2, X2 is predominantly a wino, and x5 
is a bino-dominated state. The decay of x\ Xi an d a Z° is kinematically disallowed. 
The branching ratio of the three-body decay is shown in Fig. 0] as a function of tan /? for 
the singlet well as for the representation 75. We have calculated the branching 

ratio using the program SDECAY (27 1. In this figure the initial value of tan/5 is 4.5, since 
for a lower value of tan/? the light Higgs mass is less than 114.4 GeV, which is the 
LEP lower limit We see from the figure that for higher values of tan (3 this branching 
ratio decreases since the branching ratio x% ~^ Xi T+T ~ increases with tan/? due to a larger 
Yukawa coupling. 

For the representation 75, X2 ls wino-dominated and Xi is bino-dominated as in the singlet 
case. However, the mass difference between the X2 an d the Xi is much smaller compared 
to the singlet case. As we see from Fig. in the low tan/? region the branching ratio 
for these two different representations are very close though the branching ratio for the 75 
representation is slightly larger. This is due to the fact that BR(x2 —* XiQQ) is slightly larger 
in the singlet case as compared to the case of 75 dimensional representation. The leptonic 
branching ratio is then almost equally distributed among the available channels. However, 
for large tan/? the branching ratio in the Xi r+r ~ channel is larger for the 75 case than 
for the singlet case. For large tan/? this makes the branching ratio in the xV- + l~ channel 
smaller for the case of 75 dimensional representation. We also note that in the case of 75 
dimensional representation the partial decay width of X2 Xi u ^ is larger than the partial 
decay width of X2 ~^ Xil + l~ m the large tan /? region. On the other hand, in the singlet case 
the partial decay width of X2 ~^ x\ v ^ is always smaller than that of X2 ~^ Xil + l~- 

For the set of parameters that we have discussed and in the case of the representation 
200, the spectrum is such that all the left and right handed sleptons are lighter than X2 an d 
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FIG. 5: Branching ratio as a function of tan 0, in the case of < m^o , for the representation 
200. All other parameters are same as in Fig. 

are produced on mass-shell. Although, Ir proceeds with 100% branching ratio to Xi + h m 
the case of li one should multiply by the appropriate branching fraction. This is shown in 
Fig. 03 One can see that the dependence on tan f3 is not significant in this case. The results 
for the representation 24 are not shown here due to the fact that it results in the lightest 
neutralino mass below the current experimental lower limit. 



1. The case of 200 

In this subsection we will consider the representation 200 of SU(5), as it arises in (jlj), in 
some detail. The ratio of the U(l) and the SU(2) gaugino masses is approximately given 
by (at the one- loop level 1 ) 

|Mi| : \M 2 \ = 2.5 : 1. (17) 

This resembles very much the scenario of anomaly mediated supersymmetry breaking for 
values of /x larger than Let us now highlight two important characteristics of this 

1 Using the two-loop renormalization group equations (RGEs) the ratio is approximately 2.6:1 for a wide 
range of parameter choices. Similarly, for other representations the change in this ratio with the use of 
two-loop RGEs does not change our conclusions. 
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representation. 

• Xi an d Xi are almost exclusively winos, and they are nearly degenerate in mass. 

• X2 1S predominantly a bino for \fi\ > Mi. 

Consider the decay X2 ~ > xV- + l~ f° r the 200 dimensional representation. We again choose 
the scalar masses in such a way that m^, m fl < m^o < m^, mp. The reason for such a 
choice is that the two-body decay X2 ~~ > W is allowed. Although the decay of Ir into I + Xi 
is highly suppressed due to the very small bino component in Xii Ir will decay eventually 
in this mode with a one hundred percent branching ratio. Of course, one should be careful 
to consider the possibility of a displaced vertex in the decay of Ir. The BR(x2 — > Xi^ + ^ ) 
calculated in this manner depends very strongly on /i (increases with increasing //) since 
as fi increases X2 becomes more and more bino like, and thus the partial decay widths of 
X2 Xi W T and X2 ~^ xV 1 are suppressed and the partial decay width of X2 W is 
enhanced. This makes the branching ratio into IrI mode larger for large values of fi. This 
is shown in Fig® The branching ratio in the channel x9> ~> Xi% is always very small for all 
values of tan (3. 

We note that in the case of the 200 dimensional representation we use the constraint m -± 
> 88 GeV applicable for nearly mass degenerate lighter chargino and the lightest neutralino 

H- 

Let us now discuss the BR(%2 — * xV- + l~) f° r this set of parameters for the representa- 
tions 1 and 75. We do not compare the case for the representation 24 here since for the 
parameter choice of this figure the 24 dimensional representation always produces a lightest 
neutralino with mass below the current experimental limit For this set of parameters 
the representations 1 and 75 give similar kind of spectrum so that no two-body decays of 
X2 are allowed. In Fig. [7| we show the branching ratio for these two representations as a 
function of tan j3 and for a value of fi = 500 GeV. 

2. The case of 24 

In this subsection we will consider the case of 24 dimensional representation where \Mi\ w 
0.166|M2|. We look for a set of parameters such that the mass of the lightest neutralino is not 
below the current experimental lower limit as was the case in the previous sub-sections. For 
the present study we again consider the mass spectrum m^o < m~ lR < m^o. We have taken 
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FIG. 6: Branching ratio as a function of tan /3, in the case of < m^o and for the representation 
200 and for three different values of fi. 

M.2 = 750 GeV and \i = -200 GeV. For this choice of the parameters the lightest neutralino 
Xi is mostly a bino with some higgsino admixture whereas the second lightest neutralino 
is higgsino dominated with very small wino and bino components. The soft masses for the 
left-sleptons are assumed to be 300 GeV, whereas those of the right-sleptons are taken to be 
150 GeV. Other parameters such as squark masses and the trilinear scalar coupling A t are 
same as before. Once again we have taken the values of these parameters at the electroweak 
scale. With this set of parameters the following two-body decay channels are dominant : 
X2 ~* an d X2 ~* ^i r - We have plotted the branching ratio of X2 ~^ Xil + l~ as a function 
of tan j3 in Fig. |H1 We see from this Fig. that for some values in the low tan /3 region this 
branching ratio can be as large as 65%. For large values of tan/?, the BR(x°, ~~ * Xi T+T ~) 
dominates. Also, for this choice of parameters in the 1, 75 and 200 representations in (0J) 
we always have a stau LSP which we do not consider in our R-parity conserving scenario. 
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B. Decay of heavy Higgs bosons into a pair of neutralinos : H°,A° — ► X2X2 

In this subsection we will study the branching ratios of the heavy Higgs bosons H° and 
A into a pair of second lightest neutralinos. We have used the package HDECAY to 
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FIG. 7: Branching ratio as a function of tan (5, in the case of mj > m^o and for the representations 
1 and 75 and for /i = 500 GeV. All other parameters are the same as in Fig. HO 

calculate the branching ratios. The decay widths and the branching ratios depend on the 
ratio of Mi and M2 along with other MSSM parameters. We have calculated the branching 
ratio of H°,A° — > X2X2 f° r different SU(5) representations that arise in the product @. 
The coupling of the heavy Higgs boson H° with a pair of neutralinos is given by j^ll 1^2 1 : 

H°XiX°- -ig(A L P L + A R P R ), (18) 

where Pl = — 75) and P R = 5(1+75) are the usual projection operators. The coefficients 
of Pl and P R are given by 



A-l = Qj* cos a - S"* sin a, A R = Q'[- cos a - S% sin a, 



(19) 



where 



Qij — 2 [^3(^2 — ZjxtanOw) + Zj 3 (Z i2 — Zntan9w)}^i, 
S'ij = 2"[^ 4 (^? 2 — ZjitanOw) + ^4(^2 — Zntan6w)]u- 



(20) 



Here Z is the neutralino mixing matrix in the basis {—iB, —iW, Hi, H2), and e« is the sign of 
the i th neutralino mass eigenvalue. Furthermore, sin a and cos a are the usual Higgs mixing 
angles. 
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FIG. 8: Branching ratio as a function of tan /3, in the case of < m^o and for the representation 
24. The choice of other parameters is described in the text. 



Similarly 



given by 



31 



the coupling of the pseudoscalar Higgs boson A to a pair of neutralinos is 



3J : 



A°x°Xj- -g(B L P L -B R P R 



(21) 



where the coefficients of Pl and Pr are given by 

B L = Q]\ sin (3 - S'!* cos /3, B R = Q% sin (3 - cos (3. (22) 

As an example, in Fig. El we have shown the dependence of branching ratio BR,(H° — > 
X2X2) on m A for a particular choice of MSSM parameters. This point in the parameter space 
is the same as in Fig. 0with the choice of tan/? = 6.5. This way we can directly compare 
the branching ratios in the representations 1, 75 and 200. In Fig. we have plotted the 
branching ratio BR(A° — > X2X2) as a function of tua for the same choice of parameters, and 
for the same SU(5) representations. We can see that for this choice of the parameter set 
and for < 350 GeV, the branching ratio of the decay of A is larger than that of the 
decay of the heavy Higgs scalar H° for the representations 1 and 75. This is due to the fact 
that for H° the total decay width is larger due to the increase in the number of available 
channels to the SM particles, which leads to a smaller branching ratio to sparticles. In the 
case of 200 dimensional representation the threshold opens up for heavier rriA, and once 
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FIG. 9: The branching ratio of H° —>■ X2X2 as a function of itla for three different SU(5) rep- 
resentations in @. Here tan/3 is taken to be 6.5, and other MSSM parameters are same as in 
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FIG. 10: The branching ratio of ^4° — > X2X2 as a function of for three different SU(5) 
representations. Here, tan/3 is taken to be 6.5, and other MSSM parameters are same as in Fig. EH 

again the branching ratio of A is larger than that of the H°. As we have discussed earlier, 
the representation 24 produces a very light lightest neutralino for this choice of parameters, 
and is not further discussed here. 
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FIG. 11: Contours of a(pp — > H°,A° — > X2X2 — > 4/ + X) in fb, where I = e ± or n± and 
X represents invisible final state particles. This is the case for the singlet representation. Other 
MSSM parameters are the same as in Fig. |§J ^/s = 14 TeV. 



C. Signal cross section 



We now consider signal cross section and the total event rate in the four lepton channel 
at the LHC with ^/s = 14 TeV for two different representation, the singlet and 75. We 
show the contours of constant cross section in the (uia, tan/5) plane for a representative 
set of MSSM parameters at the electroweak scale. As in Fig. we have taken \x = 500 
GeV, M 2 = 120 GeV, all left slepton masses to be 300 GeV and all the right slepton masses 
to be 150 GeV. All squark masses are taken to be 1 TeV. The top mass is m t = 178 GeV 
and the bottom mass nib is 4.25 GeV. The production cross section gg — > H°/A° has been 
calculated in the next-to-leading order using the package HIGLU 1 331 , which is based on the 
calculations in Ref.p4J. For the gluon distribution function we have used the distribution 
given in j^. We note that for low values of tan/5 this channel dominates the production 



cross section. We have a 
at the leading order 



so considered the inclusive associated production qq, gg — > bbH°/A° 
(which is essentially the leading order subprocess bb — > H°/A°). 
The factorization and the renormalization scale are chosen to be //p = /jr = (m^/ J 4 + 2m; ) )/2. 
The process gg — > bbH°/A° as well as the process gg — ► H°/A° is enhanced for larger values 
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FIG. 12: Contours of a(pp — ► H°, A — ► X2X2 — > <±l + X) in fb, where I = e ± or // ± and X 
represents invisible final state particles. This is the case for the representation 75. Other MSSM 
parameters are the same as in Fig. El yfs = 14 TeV. 

of tan f3 due to the large coupling of Higgs bosons to bb. However, the process gg — > bbH°/A 
dominates the production process for the values of tua that we have considered here {uia ~ 
200 GeV). In order to get a quantitative idea of these individual contributions to the Higgs 
boson production let us give an example here. If we take tan f3 — 20 and = 230 GeV 
then gg —>■ H production cross section is 4.08 pb and Hbb production cross section is 31.4 pb 
at LHC energy. Next, we multiply these Higgs production cross sections by the appropriate 
branching ratios B.K.(H°/A° — > X2X2) an d B.R.(%2 "~ ¥ Xil + l~) discussed in the previous 
subsections in order to get the four lepton signal at the LHC. In Fig. El we have shown the 
contours of constant cross section of the 41 signal for the singlet representation arising in 
(HJ). We see that for tan/? up to ~ 10, and ~ 250-350 GeV the total Al cross section can 
reach up to 10 fb, which corresponds to 1000 signal events (without any cuts) for integrated 
luminosity of 100 fb -1 . In Fig. El the contours of constant cross section are shown for the 
representation 75 for the same choice of parameters. It is evident from this figure that a 
smaller region in the (m^, tan/5) plane can be probed in this case with the same number 
of events. However, it shows different possibilities for these two representations. It is also 
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evident from these two figures that the four lepton signal is very small for large values of 
tan/? due to the suppression of the branching ratio B.R.(if° /A — > x^X®)- ^ n the case of 
the representation 200, this four lepton signal is available only for higher values of (> 
600 GeV) as can be seen from Fig. |U] and Fig. El for identical choices of other parameters. 
However, the total cross section is less than 1 fb for a large region in the (m^-tan/3) plane 
and we do not show any separate plot for that. The total cross section in this four lepton 
channel can be similarly studied for the representation 24 for some different set of parameters 
which we do not pursue here. 

We now briefly discuss the possible backgrounds to this four lepton signal. There can 
be two types of backgrounds, namely the Standard Model processes leading to this type of 
signal, and SUSY processes. The main SM background comes from Z°Z° and ti production. 
As discussed in Ref. jl5|, the ti background can be eliminated to a large extent by requiring 
four isolated leptons with pf > 10 GeV. Demanding a missing transverse energy of 20 GeV 
and an explicit Z° veto can reduce the background from Z°Z°. The background from SUSY 
processes can come from squark/gluino production or sneutrino pair production. The events 
coming from squark/gluino production can be eliminated by requiring soft jets with Et < 
100 GeV and E™ 1SS < 130 GeV. The background coming from sneutrino pair production 
is more difficult to handle. However, it could possibly be distinguished due to the fact 
that it has larger E™ 1SS and larger pf compared to the signal. After using all these cuts the 
percentage of four lepton signal surviving is approximately about 60%. A detailed simulation 
of the signal and background events is beyond the scope of the present work. We hope to 
come back to this issue in a future work. 

IV. HIGGS PRODUCTION IN THE CASCADE q,g -> x° -> ^Xi 

If squarks and gluinos are light enough to be produced (pp — > qq' , gg, gq), then their 
production cross section will be large at a hadron collider. Thus, the decay chain 

q, xl + X - x°ih(H°, A°)+X^ x\bb + X (23) 

will be an important source to look for Higgs bosons at LHC in the final state bbbb + X. This 



chain has been considered at LHC with universal gaugino masses in 



12. 161. It was found 



that for suitable values of the parameters, the signal for all of the neutral Higgs bosons was 
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clearly above the background. Discovery potential for A and H° extended to 200 GeV, 
independent of the values of tan (3. 

Here we will consider the decay chain in (|23[) for nonuniversal gaugino masses, and study 
the changes that occur from the case of universal gaugino masses. We will only consider the 
cases without top squarks, which is enough to illustrate the differences that arise when the 
gaugino masses are nonuniversal. 

In proton collisions, a squark pair, squark-antisquark, squark-gluino, or a gluino pair can 
be produced. We have used PROSPINO [37] to calculate the production cross sections of 
these modes (for y/s = 14 TeV). Here we study the part of the parameter space for which 
rrig > TTiq. Then every gluino decays to a quark and the corresponding squark. If the mass 
difference of the (top) squark and gluino is large enough then gluino may also decay to a 
top-stop pair. Because we are considering only the five lightest quark flavors, we have to 
remove top-stop contribution from the results. 

In order to compare with the calculation in the universal (singlet) case in 12], we use 

n n 

the average branching ratios for particles, as used in [12J], and which are defined in [38] . 
Thus we sum over all the decay widths of squarks decaying into a quark and a neutralino 
divided by the total decay width of squarks decaying into any neutralino and quark. The five 
squarks (excluding the stop) are considered to be equal in the sense that the bottom Yukawa 
coupling effect is neglected. The decay branching ratios are calculated using SDECAY [27 1 
for squarks and neutralinos, and HDECAY J^iJ for the final decay of Higgs bosons to bottom 
quarks. 

In Fig. (a) we have plotted the cross section of the decay chain ()23)1 as a function of the 
gluino mass for the singlet representation. We have used parameters tan (3 = 10, = 200 
GeV, fj, = +500 GeV, nig = 600 GeV and mj = 350 GeV as low scale input values. 

In the case of singlet representation, only the decay through the light Higgs boson h° 
is kinematically possible. The mass difference of the two lightest neutralinos is too small 
to produce heavier Higgs bosons H° and A . We can see a sharp rise in the cross section 
where the light Higgs channel opens up. This is due to an increase in the X2 an d Xi mass 
difference as a function of the gluino mass. The production cross section of squarks decreases 
as gluino mass increases. This is of course independent of the representations arising in the 
product (J1J). 

In Fig. (b) we have plotted the corresponding cross section for the 24 dimensional 
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FIG. 13: Cross section pp — > bbbb +X at LHC through the decay chain q, g — > x[> ~* x\h(H° , A ) 
x5&6 + -X" for (a) the representation 1 via /i, and (b) the representation 24 via h, H° and A . 



representation. Because of the changed relations between the gaugino mass parameters, 
the composition and masses of the neutralinos are different from the universal case. Now 
all the Higgs channels are available. We see that the CP-even neutral Higgs if°-channel 
gives the largest cross section. As gluino mass increases, the decay branching ratio of 
X2 —* h, H°, A + Xi decreases. 

In the 75 dimensional representation the mass difference of the two lightest neutralinos 
are generally too small in order for %° to decay into Higgs bosons, thereby making the decay 
chain, Eq. (|2"3*)l . irrelevant. 

For the 200 dimensional representation, the mass difference of the two lightest neutralinos 
depends on the squark masses. Requiring the lightest neutralino to be the LSP, and for 
experimentally viable Higgs boson mass, the mass difference of the two lightest neutralinos 
is relatively small, and the total cross section resulting from the decay chain (|23|) remains 
below the detection level. 

V. CONCLUSIONS 

We have studied the consequences of gaugino mass nonuniversality as it arises in a su- 
persymmetric grand unified theory for neutralino masses and mass relations, as well as for 
particular Higgs production and decay processes. 

We found that the upper bounds of neutralino masses and the mass sum rules depend 
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significantly on the representation. Similarly the studied decay possibilities of Higgs bosons 
depend on the representations. The decay of the second lightest neutralino to two leptons 
and the lightest neutralino very much depends on the mass difference between the lightest 
and second lightest neutralino which in turn depends on the representations. This is also 
true for the production of Higgs bosons in the decay of the second lightest neutralino. 

From our considerations it seems clear that depending on the region of the parameter 
space, the Higgs decay h(H°, A ) — > X2X2 ma y be observable for the gauginos emerging 
in any of the representations 1, 24, 75, or 200 of SU(5). However, the region in which 
X2 — > 21 + X is large, and it is possible for Higgs bosons to decay to the second lightest 
neutralinos, is rather limited in any of these models. Thus, the relevant regions of the 
parameter space do not necessarily overlap. In Section 3 we compared singlet and 75 using 
the same set of parameters. For 200, the total cross section is less than 1 fb for the same 
choice of parameters when kinematically available. For the representation 24, we did not 
find a region where a comparison could have been made. 

Interestingly, for the production of the Higgs bosons via the decay chain including 
h(H , A°)xi, in addition to the singlet, we found relevant region of the parameter space 
only for the representation 24. Furthermore, in this region the signal cross section for both 
neutral heavy Higgs bosons is reasonably large for not very heavy gluinos. It should be 
noted that in these two cases the signatures are clearly different. In the representation 24 
the cross section is largest at the lighter values of the gluino mass for the heavy Higgses. For 
the lightest Higgs the production channel is open for all of the discussed gluino masses in 
the representation 24 as opposed to the case of 1 representation. Also the fact that all the 
neutral Higgs channels are open in the 24 case distinguishes it from the singlet case, where 
only the light Higgs channel is available. 

Finally, we note that it is p ossible to find similar signatures from the scenarios with 
non- universal Higgs masses [39j ]. 
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